Kramers-Moyall cumulant expansion for the probability distribution of parallel 
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A general equation for the probability distribution of parallel transporters on the gauge group 
manifold is derived using the cumulant expansion theorem. This equation is shown to have a 
general form known as the Kramers-Moyall cumulant expansion in the theory of random walks, 
the coefficients of the expansion being directly related to nonperturbative cumulants of the shifted 
curvature tensor. In the limit of a gaussian-dominated QCD vacuum the obtained equation reduces 
to the well-known heat kernel equation on the group manifold. 
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I. INTRODUCTION 

One of the most popular test objects used to detect 
quark confinement is a Wilson loop calculated in some 
representation of the gauge group. Measurements of 
the Wilson loops allow one to investigate the interac- 
tion potential between static quarks in different repre- 
sentations of the gauge group. For instance, Wilson 
area law corresponds to QCD string with constant ten- 
sion 0. Most numerical simulations indicate that Wil- 
son loops indeed obey the Wilson area law at interme- 
diate distances ~ 0.2 ... 1 fm 0, 01- However, at large 
distances QCD string may be torn apart by a quark- 
antiquark pair born from the vacuum, or the charge of 
quarks may be screened by dynamical charges, and as a 
result Wilson area law is violated. Another distinct fea- 
ture of QCD vacuum is the Casimir scaling phenomenon, 
which was observed in lattice simulation 0, Q ■ Violation 
of Casimir scaling is also related to screening of static 
sources. Screening effects in SU{2) Yang-Mills theory 
were detected in numerical simulations 0,0- A proper 
theoretical description of screening in QCD vacuum is 
still absent, although many beautiful phenomenological 
models were proposed (for a review see, for example, 
00). 

An interesting phenomenological description of screen- 
ing effects, which is based on the mathematical anal- 
ogy between loop dynamics and the motion of a random 
walker on the g auge gro up manifold, was proposed re- 
cently in H 0, lift ITl| . The theory of random walks 
on Lie groups is an interesting mathematical subject 
[l2l fl3| , and one can expect that physical applications 
of the results obtained therein may provide important 
insights into the physics of confinement. For example, 
an important results concerning the connection between 
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Chcrn-Simons theory, B-model branes on noncommuta- 
tive two-sphere and random walks in the fundamental 
Weyl chamber of Lie algebras were obtained in . In 
the work Q it was observed that Casimir scaling corre- 
sponds to a free diffusion on the group manifold, which 
is described by the heat kernel equation. This could be 
expected, since it is known from lattice gauge theory 
that infrared behavior of Yang-Mills theories in higher 
dimensions can be described by two-dimensional Yang- 
Mills theory (dimensional reduction scenario, 14]), and 
the partition function of two-dimensional Yang-Mills the- 
ory obeys the heat kernel equation In order to de- 
scribe screening effects one should somehow modify the 
simplest heat kernel equation. A straightforward modifi- 
cation is to add a drift term to the equation, which was 
interpreted as the Fokker-Planck equation on the group 
manifold in 0, [j| . Such a modification led to many non- 
trivial physical consequences, such as the emergence of 
confining fc-strings Zn symmetry of the group cen- 
ter was found to play a crucial role in transition from 
Casimir scaling to screening regime Q , which has a 
nice interpretation in terms of central vortices. A non- 
trivial but solvable modification of QCD2 matrix model 
which describes screening effects was discovered in 
In our work |10| we have pointed out another possibility 
to generalize the heat kernel equation, namely to include 
the terms with higher-order derivatives with respect to 
group coordinates. An advantage of such a generaliza- 
tion is that there is no need to introduce some external 
force in the Fokker-Planck equation " by hands" , as it is 
done in the works H,0 ■ Thus geometric properties of the 
group manifold (homogeneity and absolute parallelism) 
are preserved and the symmetry of the equation under 
group transformations is explicit. 

The aim of this paper is to provide a theoretical ground 
for the phenomenological analysis of 0, 0, and to de- 
rive the most general differential equation for the prob- 
ability distribution of parallel transporters on the group 
manifold from the first principles. The classical loop 
equation [Tfjl ] is first rewritten as a flow equation on the 
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group manifold. Then we apply the cumulant expansion 
theorem to this flow equation and find that the proba- 
bility distribution obeys the equation which can be ob- 
tained from the heat kernel equation by adding higher- 
order differential operators. This equation generalizes 
the Fokker-Planck equation used in [a 13 and is known as 
the Kramers- Moyall cumulant expansion in the theory of 
random walks |17| . The coefficients in the obtained equa- 
tion can be expressed in terms of nonperturbative cumu- 
lants of the shifted curvature tensor, which are the basic 
objects in the method of field correlators 0,0,0]. In 
the limit of a gaussian-dominated QCD vacuum, where 
only the second-order cumulant survives, this equation 
reduces to the heat kernel equation, in accordance with 
the results of Hcj . We suppose that the proposed 
generalization of the heat kernel equation may be help- 
ful in finding effective action for QCD string which takes 
into account screening effects. Kramers-Moyall cumu- 
lant expansion also emerges naturally as a continuum 
approximation of discrete random walks, which can be 
interpreted as a continuum limit of a lattice gauge the- 
ory. 

The paper is organized as follows: in this section we 
introduce the loop variables and the classical loop equa- 
tions 0]. Since we extensively use geometric construc- 
tions on the group manifold throughout the paper, in the 
next section we make a brief review of geometric prop- 
erties of compact simple groups and introduce the group 
laplacian. In the section ITTT1 the probability distribution 
of parallel transporters on the group manifold is intro- 
duced and its classical limit is considered. Finally, in the 
section llVl we apply the cumulant expansion theorem to 
obtain the main result of this paper, namely the differ- 
ential equation for the probability distribution of parallel 
transporters on the group manifold. 

We consider gauge fields A^(x) in some representa- 
tion of the Lie algebra of the gauge group G, which is 
supposed to be compact and simple. To any smooth 
closed path 7 one can attribute parallel transport op- 
erator which can be written as a path-ordered exponent: 



# (7) = 73 exp 



dx^A, 



(1) 



Parallel transport operator is the element of some repre- 
sentation of the group G, therefore in classical theory to 
each loop 7 corresponds some element 50(7) of the gauge 
group G. In order to calculate the path-ordered expo- 
nent in one should choose some initial point xq on 
the loop 7. Parallel transport operator transforms under 
gauge transformations in the following way: 



# ( 7 )^f (/(*„))* ( 7 ) f (rHxo)) 

Soil) -»• f(xo)go(l)f' 1 (x ) 



(2) 



be identified, and points in the physical phase space are 
actually the group classes. 

We consider a one-parameter family of loops 7 (r), as- 
suming that at fixed r the loop 7(7") is parametrized with 
some variable a, i.e. 7 (r) : a — > x^(a,r), x M (0) — 
x tJ, (2n). Group element which corresponds to the loop 
7(7") in classical theory will be denoted as go(r) = 
.9o (7( T ))- ln order to impose boundary conditions in sub- 
sequent equations we suppose that r = corresponds to 
a loop of zero area, for which g (7(0)) = 1. We will also 
assume that the worldsheet swept by the loops 7(7-) is 
simply connected. This assumption allows us to intro- 
duce the shifted curvature tensor F^ v on this worldsheet, 
as is done in the proof of the non-Abelian Stokes theorem 
MM: 



Fp V {x) = U(x ,x,x)F^(x)U(x,x ,x) 



(3) 



where F^ v = d^A^ — i[A^,A u ] is the curvature tensor 

/ x \ 

of Yang-Mills field, U(x,y,x) = "Pexp i J dx^A^ is 

\ y;x ) 

the non-Abelian phase factor along the path x, xq is the 
reference point on the surface and x is some path which 
connects the reference point x$ and the point x. Shifted 
curvature tensor is the curvature tensor parallel trans- 
ported from the point x to the reference point xq along 
the path %. 

We start with the Polyakov-Migdal variational equa- 
tion for loop variables in classical Yang-Mills field fjlq] . 
chapter 7): 



5*$ fa) ? +/ , ~ , , dx v 



(4) 



Thus gauge orbits for loop variables are group classes, i.e. 
the points which belong to the same group class should 



For our continuous set of loops the equation reduces 
to a linear differential equation: 

= i J dad a x u (a,T)d T x fJ -(a,T)F^ (x(a,r)) 

7(t) 

where *(r) = * (7(7-)) and d a = -§^, d T = Jp Boundary 
condition for this equation is ^(0) = /, because we have 
assumed that go (0) = 1. 

II. GEOMETRY OF THE GROUP MANIFOLD 

In order to proceed with the analysis of the equation 
(JSJ it is necessary to introduce some geometric construc- 
tions on the group manifold. In this section we make a 
brief review of geometric properties of simple compact 
groups. In order to relate our notations to those of an- 
other authors 0, we also introduce here the group 
laplacian. 

Group generators in the fundamental representation 
T a , — T a are chosen to satisfy the normalization 
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Tr T a Ti, = and obey the commutation relations 
[T a ,Tf,] = iC^ h T c . Suppose that group elements are 
parametrized by some coordinates g a (for the sake of 
brevity we will use the same letters to denote the ele- 
ments of the group and their coordinates). In this paper 
we will work in several different spaces: the space-time 
where our Yang-Mills fields live, the group manifold and 
the Lie algebra (tangent space to the group manifold at 
the group identity). From now on we will use latin in- 
dices a, b, c, . . . to denote the elements of the Lie algebra, 
the first greek indices a, f3, . . . for tensors on the group 
manifold, and the indices from the middle of the greek al- 
phabet /i, v, . . . for tensor fields in space-time. Wc define 
the vector fields L" (g) and R" which are induced by the 
left and right action of the group generators respectively 



infinitesimal version of this constraint is: 



{l + ie a T a )U{g a ) = U(g a + LZ{g)e a ) 
U(g a )(l + ie a f a ) = U(g a + R a a (g)e a ) 



(6) 



where U(g a ) are the matrices of the fundamental rep- 
resentation and e a are arbitrary infinitely small param- 
eters. Let us also introduce the derivatives V a and V a 
(not to be confused with Yang-Mills covariant derivatives 
in space-time) along these vector fields: 



(7) 



These derivatives obey the following commutation rela- 
tions: 



[V a ,V 6 ] = c a c h v c 



v a ,v t 



= -c c ab v c 



V„,Vt 







(8) 



Thus the derivatives V a or V a acting on the functions 
on the group manifold build an infinite-dimensional re- 
ducible representation of the Lie algebra of the group G. 
Corresponding representation of the group G is the reg- 
ular representation, which contains all finite-dimensional 
irreducible unitary representations. Quadratic Casimir 
operator of this representation is |24| : 



A = V a V Q = V a V a 



(9) 



Eigenvalues of —A are the eigenvalues of the quadratic 
Casimir operators in irreducible unitary representations. 
The minus sign arises here because the derivatives V a 
are anti-hermitean, unlike the hermitean generators T a . 

While the derivatives V a and V a generate left and 
right group multiplications respectively, their difference 
generates shifts within group classes: 

V a L>(.g) - if a U{g), V a U{g) = iU(g)f a 

V a - V a ) U(g) = i [f a , U(g)] (10) 

In this paper we will work in a Hilbert space of func- 
tions on group classes. Such functions obey the con- 
straint 4>{g) = 0(/g/ -1 ) f° r an y group element /. An 



(v„ - v.) m = o 



(ii) 



Thus the Hilbert space of functions on group classes is 
the subspace of the full Hilbert space of functions on the 
group manifold, spanned on the solutions of the linear 
equation 

An invariant metric on the group manifold is fixed up 
to a constant by requiring infinitesimal distances to be 
invariant under group multiplication |22j j. This means 
that the vectors L a a are Killing vectors, which immedi- 
ately gives for the group metric (Killing form) h a p\ 



h°P(g) = L«M(g), h^hp 1 = 81 



(12) 



An invariant measure on the group manifold (the Haar 
measure) is dfi(g) — yj det h a p(g) L] dg 1 ' . 

7 

Group structure of the manifold allows also to define 
the parallel transport and the connection on it. Paral- 
lel transport is defined as follows: consider two infinitely 
close points with coordinates g a and g a + dg a . Sup- 
pose that under left group multiplication these points 
transform into g a and g a + dg a . Then dg a is said to 
be the vector dq a parallel transported from the point 
g to the point g |22| . The corresponding connection is 
= LJ-j^Lp a , where indices are lowered and raised 
with the group metrics h a p. One uses this connection 
to define a covariant derivative V a on the group mani- 
fold. Group manifold is a space of absolute parallelism, 
i.e. the space with zero curvature but nonzero torsion: 
[V a , Vp] = SlpVj. We use the covariant derivatives V a 
to redefine A and to extend its action on arbitrary tensor 
fields on the group manifold: 



A = h af} V n V 







(13) 



For scalar functions the definitions 10 and (|13fl are equiv- 
alent, as in this case A is the Beltrami-Laplace operator 
for the group metrics h a p: 

A = (det Kp)- 1 ' 2 -^ ((det Kpfl^-^j (14) 

Definitions (|13fl and l|14fl are different in the case of ar- 
bitrary geometry, but in the particular case of manifolds 
with group structure they coincide. 



III. PROBABILITY DISTRIBUTION OF 
PARALLEL TRANSPORTERS ON THE GROUP 
MANIFOLD AND ITS CLASSICAL LIMIT 

In quantum theory one can not establish any deter- 
ministic correspondence between loops in space-time and 
points on the group manifold. Instead one can consider 
the probability distribution p(g,r) which will determine 
the probability dP(g, r) = p(g, T)d/J,(g) for the group el- 
ement g(j) associated with the loop 7(t) to be within 
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an infinitesimal volume dfi(g) on the group manifold. As 
the physical phase space contains only group classes, the 
points which belong to the same group class should be 
identified, which means that the physical probability dis- 
tribution is actually defined on group classes. Thus we 
require the function p(g, r) to be the function on group 
classes, i.e. to obey the constraint (|1H . As will be 
shown further, the constraint is consistent with the 
final equation for p(g,r). A full orthonormal basis in 
the Hilbert space of functions on group classes is built 
by the characters of irreducible unitary representations 
Xk(g), where the index k labels irreducible unitary rep- 
resentations of the group G, including the trivial one. An 
identity operator on this space is the delta-function on 
group classes: 



*(/,<?) =X>*(/)x*(ff) 



(15) 



where we have taken into account that 5 (g, f) = 5 (/, g) 
and hence -J^ 5 (g, /) = — JL? S (g, /). The vector field 



dg d ^ can be found from (JSJ: 



dgS(r) 
dr 



L a a (go{r)) J dad^d^F^ix^r)) (20) 

7(r) 



It is convenient to introduce the linear differential oper- 
ator C(t) which acts on functions on the group manifold 

as the derivative along the vector field — dg ° — • 



At 



C(t) = - j dad a x v d T x»F* v {x(a, r)) V a (21) 

7 (t) 

Now the equation 119(1 can be rewritten as: 



By definition, for any function on group classes 4>(g) one 
has: 



d^{g)mHLg) = <t>{f) 



As usual in quantum theory, we introduce the proba- 
bility distribution p(g, r) as the expectation value of the 
delta-function on the physical phase space 0, El El : 



p{g,r) = ( S (g,go(r)) ) 



(17) 



where by (...) we denote the vacuum expectation 
value. As (1) = 1, the probability distribution H17|) 
is automatically normalized to unity: J d/i(g)p(g, r) = 

(Jdn(g)S(g,g (T))) = l. 

G 

As S(g,g (T)) in lfT7|l is the delta- function on group 
classes (|15p. the function p(g,r) is a gauge-invariant 
quantity. A simple calculation shows that indeed the 
probability distribution p(g, t) can be expressed in terms 
of Wilson loops |icj : 

P(g, r) = {5 (g, g (r)) ) = ( ^ Xk(g)Xk (so(r)) ) = 



k 



^ Xk(g)(Xk (ffo(r)) } = Xk(g)d k W k ( 7 (r))(18) 

k k 

where W k (7) is the Wilson loop in the fc-th representation 
of the gauge group and d k is the dimensionality of the k- 
th representation. 

We would like to find a differential equation satisfied 
by p(g,r). To this end we consider the classical limit 
of p(g,r), which is given by the delta- function at the 
classical solution: po(g,r) = 5 {g,go (r)). Differentiating 
Po(g,T) over r gives: 

d T S(g,go(r)) = -^S(g,g {r)) MM = 

~£ stoM) ^ (19) 



d TPa (g,r) = -^lJL po (g :T ) = C(r)p (g.T) (22) 



(lg) which is nothing but the flow equation along the vec- 

Initial condition for this equation is 



tor field 



d 9o(T) 
dr 



Po(9,0) = as <? (0) - 1. 



IV. CUMULANT EXPANSION 

The probability distribution p(g, r) is the vacuum ex- 
pectation value of the function po(g,r) = 8(g,go(T)) as 
the functional of the gauge field A fJi (x), which enters 
the equation for Po(g,T~) through the differential oper- 
ator C{t). Our aim is to obtain an equation for p(g,r), 
which can be done without direct calculation of p{g, r) 
by using the cumulant expansion theorem (also known 
as the Van-Kampen expansion) pM l2r| . Cumulant ex- 
pansion was originally devised to solve linear differential 
equations with random coefficients. Together with the 
non-Abelian Stokes theorem it was successfully applied 
in nonperturbative QCD to find vacuum expectation val- 
ues of Wilson loops [|| 0, . An approach based on 
the cumulant expansion was further developed into the 
method of field correlators 0, fl8L IT^ | . A particular case 
of the cumulant expansion is known in the the ory of ran- 
dom walks as the Kramers-Moyall expansion 17] . 

Applying the cumulant expansion theorem to the equa- 
tion (1221 yields the following equation for p(g, r) = 

(po(g,r) ): 



d T p(g, t) = Y] J ' 

■ ■ T 



fc=0 



j ...j dn...drk 

>T 1 >T 2 >...>T k 



{{C(T)C(r 1 )...C(r k )))p(g,r) 



(23) 



where double parentheses ((...)) denote the ordered cu- 
mulants [2{J , which extend the concept of cumulants 
of scalar functions onto non-commuting operators (note 
that the operators £(r) for different values of r do not 
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commute). Ordered cumulants (( ... }} are defined im- 
plicitly by the following recurrence relation |25l |26| : 

n 

(£(r)£(n)...£(r n ))=^ £ (24) 

r=0 {„..., r }[J{j 1 ...) I1 _ r )={l...») 

((£(r)£(r n ) . . . C(r ir )}) (C(r h ) . ..L(r jn _ r )) 

where summation is performed over all decompositions 
of the set {l...n} into two ordered subsets {i%...i r } 
and {ji . . . j n - r }, each of which may be empty. To start 
the recursion we define (( C(t) )) = ( £(t) }. Ordered 



cumulants of the lowest orders are |25j: 

((£(n)£(T 2 )» = (£(n)£(r 2 )) - (£(n)){/:(T 2 )) 
«£(n)£(T 2 )£(T 3 ))) = (£(T 1 )£(r 2 )£(r 3 )> - 
-(£(n)£(T 2 )><£(r3)) - (£(r 1 ))(£(r 2 )£(r 3 )) - 
-(£(r 1 )£(r 3 ))<£(r 2 )) + ^(n))^))^)) + 

+(£(r 1 ))(£(r 3 ))(£(r 2 )) (25) 

Explicit substitution of £(r) in (|23[1 gives the desired 
equation for p(g, t): 



d T p(g,r) = ^(-l) fc+1 /" dad^d^ / ... j ' dS" lUl ...dS^ k 
k=0 7(r) 

(( F;, (x(a, r)) ^ (x(<7i,n)) . . . F^ k (x(a k ,r k )) ))V„V 0l . . . V Qfc p(. 9 , r) 



(26) 



where = l/2dadrd a x^ d T x v ^ is the volume form on 
the worldsheet swept by the loops 7(7") and integration is 
performed over the whole worldsheet with the restriction 
t > n > t 2 > . . . > Tfc . Note that in the above equation 
we have used the cumulants of the shifted curvature ten- 
sor in the basis of the generators T a . These cumulants 
are defined by the following relation: 



-((F^ (zfa, n)). 



F% Vk (x(a k ,T k ))T ak )) 



■F% Vk (x(a k ,T k ))))T ai 



view of the properties of such cumulants can be found in 
For instance, the lowest-order cumulants are: 

//p« pb )) = (p a F b ) — (F a )(F b ) 

W- 1 Ml l>i * M2^2'' \ Ml "1 M2^2' \ Ml^l/ \ M2 W 2 / 



((F a F b F c )) = (F a F 

W" Ml^l M2f2 IJ.3V3II \ Ml"! / 



F c 

M2^2 M3"3 



Ml f 1 M2 ^2 M3 ^3 ' ' \ M 

-(F a F b )(F C ) - (F a )(F b F 

\ Ml"! M2f2 ' \ M3^3 ' > Ml"l ' » M2^2 M3^3 



-(F a )(F b F c ) + (F a )(F b )(F C ) + 

X- 1 Ml^l ' ' M3^3 M2^2 ' 1 \ Ml^l ' > M2"2 / \ f>S"3l 

+ (K lVl )(F b ^ 3 )(F^ 2 ) (28) 



where {{F^ lVl . . . F^ kVk )) are the nonperturbative cumu- 
lants of the shifted curvature tensor |6j, Il8l Il9| . Such cu- 
mulants are defined by the relation similar to (|24|l . where 
r-ordering is replaced by some surface ordering |20l l2ll| , 
which we will assume to be the ordering with respect to r 
and a variables. It should be noted that the shifted cur- 
vature tensor always transforms as a curvature tensor in 
the reference point xq, thus the cumulants (|27|l transform 
covariantly under gauge transformations. A detailed re- 



The equation (|26[) can be rewritten in a compact form 
as: 



drP{g, t) 



fe=0 



n 



°(T)V a V ai ...V akP (g,T) (29) 



where the coefficients rj a nfc (r) are some functions of 
t which can be expressed in terms of the cumulants of 
the shifted curvature tensor: 



n»«i-«»( T ) = (-l)W f dad a x v d T x^ ^■■■^ dS" 1 " 1 ...dS" kVk 

J T>T!>T 2 >...>T k 

7(r) 

« F« v (x(a, r)) F«V a (xfan)) . . . Fft^ (x(a k ,r k )) )) (30) 
I 

The equation (|26|) in the form l|29[) has a general form Gauge invariance imposes rather strict constraints on 
of the Kramers-Moyall cumulant expansion known in the the tensor structure of the coefficients rf °i--- a * (r). Only 
theory of random walks |l7j . colour singlets can have non-zero vacuum expectation 
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value, which means that each monomial of the form 
J7 00l " -0 " e (r)V V 01 • ■ ■ V Q|j should be a Casimir operator. 
As the derivatives V Q build a representation of the Lie 
algebra of the gauge group, this means that each such 
monomial should commute with V Q : 



v a j = o 



(31) 



Examples of such monomials are A 2 = 
S ab 5 cd X7 a V b X7 c V d ,C adj A = C a c b V a V b V c and so on. 
It is now straightforward to show that if the relation 
(|31|l holds for all differential operators on the r.h.s. of 
(|29[l . the solution of the equation (|29l) will be the func- 
tion on group classes. Indeed, for t = p(g, t) = S(g, 1) 
is the function on group classes. From the commutation 
relations JHJ it follows that the derivatives V a commute 
with arbitrary functions of the derivatives V a , and hence 



= (r)V a V ai ...V 



whenever 



<EU holds. This implies that if (V h - V h J p{g,0) = 0, 

(Vb — p{g,T~) = for all r. Thus gauge invariance 

of the cumulants (|27[1 guarantees that the solution of 
the equation (|29|) always belongs to the Hilbert space of 
functions on group classes. 

Gauge invariance requires that the "drift term" rf 
should vanish and that rf h {r) = C(r)<5 ab , where C(r) is 
some function of r which is proportional to the second- 
order cumulant of the shifted curvature tensor. In the ap- 
proximation of a gaussian-dominated QCD vacuum one 
assumes that the only non-zero cumulant is the second- 
order cumulant 0, [Hi ■ In this case the equation lj2l3|) 
reduces to the well-known heat kernel equation, in agree- 
ment with the results of H, 0, : 



d T p(g,T) = C(T)Ap(g,T) 



(32) 



Such limit of the equation (|26[) automatically leads 
to Casimir scaling (which is not surprising, since the 
gaussian dominance conjecture was devised to explain 
Casimir scaling), while possible deviations from the Wil- 
son area law can be incorporated in the function C(r). 

General form of the equation (|26|l suggests that decay 
rates of Wilson loops in different representations of the 
gauge group should be proportional to some combination 
of Casimir operators, in full accordance with the analysis 
based on the method of field correlators [|| 0, . Note 
that such conclusion is in general not true, for instance, if 
the "drift term" _n a is not equal to zero, as was proposed 
in the works 



V. CONCLUSIONS 

In this paper we have derived the differential equation 
for the probability distribution of parallel transporters 
on the gauge group. It turned out that it is necessary 
to add the terms with higher-order derivatives to the 
simplest heat kernel equation 0, U in order to re- 
produce correctly deviations from Casimir scaling and to 



describe the screening effects correctly. Gauge invariance 
requires that all the differential operators in the r.h.s. 
of the equation 12(j|) should be Casimir operators. In 
this case the tension of QCD string is in general propor- 
tional to some combination of Casimir operators, as the 
method of field correlators predicts [(J fl8L fl9|. In con- 
trast to the phenomenological approach of jaM IH3 the 
equation (12(it contains no arbitrary potential functions 
and is thus explicitly invariant under group transforma- 
tions. Moreover, the coefficients rf ai -"°* (r) are directly 
related to nonperturbative cumulants of the shifted cur- 
vature tensor and hence the equation l|26|) contains no 
phenomenological parameters at all and is exact as long 
as the cumulant expansion converges. 

The final equation (|26|) is the Kramers-Moyall cu- 
mulant expansion, which generalizes the Fokker-Planck 
equation onto arbitrary random walks |g, la Il7| . Such 
expansion arises naturally when one considers an arbi- 
trary discrete random walk and then takes the continuum 
limit ■ In lattice theories the set of loops is a discrete 
set, and the initial equation (jSJ turns into an equation 
with discrete variable t, which corresponds to a discrete 
random walk on the group manifold, according to the 
interpretation of 0, • One can note an interesting sim- 
ilarity between the continuum limit of discrete random 
walks on Lie groups and the continuum limit of lattice 
gauge theories. In particular, in both cases continuum 
limit means that physical observables do not depend on 
time step or lattice spacing. 

As long as the screening effects are small, the higher- 
order differential operators in (|29() can be treated as a 
small perturbation of the free diffusion described by the 
equation l|32|) . As it is known that the equation H32|) is 
solved by pure two-dimensional Yang-Mills theory living 
on the surface spanned on the loop 7(7-) 10], it can be 
particularly interesting to find a small perturbation of 
two-dimensional Yang-Mills theory which allows to solve 
the general equation lj2U|) at least up to several orders in 
the perturbation strength. In the case of non-zero "drift 
term" rf such modification was found in 9J. Indeed, 
the dimensional reduction scenario relates the observ- 
ables of four-dimensional and two-dimensional theories 
[l4| , which automatically leads to Casimir scaling. Thus 
possible modifications of two-dimensional theory can be 
used to extend the applicability of the dimensional re- 
duction scenario. Such modifications may be also helpful 
in finding effective actions for QCD strings which take 
screening into account. 

Possible generalization of (|26[) onto sets of loops which 
sweep multiply connected worldsheets can be also in- 
teresting. One can also implement the methods de- 
scribed in the section [H] to describe diffusion of higher- 
order differential forms. Finally, if it was possible to ob- 
tain an equation similar to l|26|) . for example, basing on 
the Schwinger-Dyson equations for quantum Yang-Mills 
fields, this will provide some information on nonpertur- 
bative cumulants of the curvature tensor. 
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